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ABSTRACT 


Presented  is  a  new  analysis  of  some  multipath  radar  data 
measured  by  Standard  Telecommunication  Laboratories  (STL)  Ltd, 
using  an  array  antenna  and  low-flying  aircraft.  The  data  have 
previously  been  analyzed  with  the  singular  value  decomposition 
(SVD)  method  by  Haykin,  Greenlay,  and  Litva.  These  authors  were 
mainly  interested  in  whether  or  not  the  SVD  method  produced  target 
elevation  estimates  consistently  close  to  an  optical  track.  The  new 
analysis,  also  using  the  SVD  method,  has  produced  additional 
results.  Firstly,  it  demonstrates  a  practical  technique  for 
selecting  target  elevation  estimates  when  the  target  range  was  less 
than  6.5  km.  Secondly,  it  provides  an  explanation,  supported  with 
simulations,  for  the  presence  of  many  false  tracks  branching  out 
from  the  target  and  image  tracks. 


RESUME 

Ce  rapport  presente  une  nouvelle  analyse  de  donn^es  radar 
concernant  le  phenomSne  de  trajets  multiples.  Ces  donnees,  obtenues 
avec  une  antenne  reseau  et  un  v6hicule  aSrien  volant  a  basse 
altitude,  ont  St4  analysees  par  Haykin,  Greenlay  et  Litva  a  I'aide 
de  la  m^thode  de  Decomposition  en  Valeurs  Singulifires  (DVS) .  Ces 
auteurs  s ' intferessaient  principalement  a  savoir  si  la  mithode  DVS 
produisait  des  estimis  de  I'felivation  de  la  cible  qui  collaient 
sans  exception  S  une  piste  optique.  La  nouvelle  technique  presentee 
ici  utilise  aussi  la  methode  DVS  mais  apporte  des  resultats 
additionnels.  En  premier  lieu,  ce  rapport  d§montre  une  technique 
pratique  pour  la  selection  des  estimis  de  1' elevation  de  la  cible 
quand  la  distance  de  celle-ci  est  infer ieure  S  6.5  Km.  En  second 
lieu,  ce  rapport  fournit  une  explication,  appuyee  par  des 
simulations,  de  la  presence  de  plusieurs  fausses  pistes  sortant  de 
la  piste  de  la  cible  et  de  celle  de  1' image. 
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EXECUTIVE  SUMMAKY 


This  report  contains  a  new  analysis  of  the  low-angle 
tracking  radar  data  measured  by  Standard  Telecommunications 
Laboratories  in  the  summer  of  1981.  The  radar  was  an  X-band  array 
antenna  located  at  the  Fraser  Gunnery  Range  near  Portsmouth,  United 
Kingdom.  The  targets  were  Hunter  and  Canberra  aircraft  flying  low 
over  the  English  Channel.  An  optical  tracker,  mounted  close  to  the 
receive  array,  provided  an  independent  measure  of  the  target 
elevation. 

The  data  have  previously  been  analyzed  by  Haykin,  Greenlay, 
and  Litva.  They  studied  whether  or  not  the  SVD  method  produced 
target  elevation  estimates  consistently  close  to  the  optical  track. 
In  the  study,  these  estimates  were  identified  as  the  candidate 
signal  elevation  estimates  closest  to  the  optical  tracks. 

The  analysis  was  very  limited  in  scope.  Besides,  it  did  not 
account  for  the  fact  that  the  SVD  method  produces  more  candidate 
signal  elevation  estimates  than  other  direction-finding  methods. 

This  new  analysis  is  concerned  with  the  development  of  a 
practical  method  for  selecting  target  elevation  estimates.  It  also 
studies  all  the  candidate  signal  elevation  estimates  generated  by 
the  SVD  method.  The  results  show  that  (a)  the  optical  track  need 
not  be  known  in  the  selection  of  the  target  elevation  estimates,  if 
the  target  range  was  less  than  6.0  km,  (b)  there  are  many  false 
tracks  branching  out  from  the  true  target  and  image  tracks,  and  (c) 
target  and  image  tracks  with  these  branches  can  be  generated  by 
adding  calibration  phase  errors  to  uncorrupted  data. 
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1  Signal  elevation  estimates  calculated  with  the  singular 
value  decomposition  method  and  the  first  STL  data  file. 
The  solid  line  is  derived  from  the  optical  tracker 
output . 

2  Signal  elevation  estimates  calculated  with  the  singular 
value  decomposition  method  and  the  second  STL  data 
file. 

3  Signal  elevation  estimates  calculated  with  the  singular 
value  decomposition  method  and  the  third  STL  data  file. 

4  Signal  elevation  estimates  calculated  with  the  singular 
value  decomposition  method  and  the  fourth  STL  data 
file. 

5  Signal  elevation  estimates  calculated  with  the  singular 
value  decomposition  method  and  the  fifth  STL  data  file. 

6  The  larger  of  the  two  signal  elevation  estimates 
derived  from  the  first  STL  data  file. 

7  Candidate  signal  elevation  estimates  derived  from 
simulated  data.  SNR  for  10  km  target  range  is  20  dB. 

8  Candidate  signal  elevation  estimates  derived  from  the 
simulated  data  used  in  Fig.  7.  Table  1  for  calibration 
phase  errors  applied. 

9  Candidate  signal  elevation  estimates  derived  from  the 
simulated  data  used  in  Fig.  7.  The  calibration  phase 
errors  are  the  negative  of  those  in  Table  1. 

10  Signal  elevation  estimates  derived  from  simulated  data. 
SNR  for  10  km  target  range  is  6  dB.  Table  1  has  been 
applied  to  the  data. 
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1 . 0  INTRODUCTION 


The  multipath  radar  data  in  this  analysis  were  obtained  by 
Pearson  and  Waddoup  [1]  of  Standard  Telecommunication  Laboratories 
(STL)  Ltd  with  a  vertically  oriented  X-band  array  antenna.  The 
targets  were  aircraft  flying  low  over  the  English  Channel.  An 
optical  tracker  mounted  close  to  the  antenna  provided  an 
independent  measure  of  the  aircraft  elevation. 

A  previous  analysis  of  these  data,  using  the 
Tufts-Kumaresan  singular  value  decomposition  (SVD)  method  [2,3], 
has  been  reported  by  Haykin,  Greenlay,  and  Litva  [4].  The  primary 
question  addressed  was  whether  or  not  the  SVD  method  consistently 
produced  target  elevation  estimates  close  to  the  true  target 
elevation.  These  authors  identified  the  candidate  signal  elevation 
estimates  closest  to  the  optical  track  as  the  target  elevation 
estimates. 

This  new  analysis  of  the  STL  data  is  concerned  with  the 
development  of  a  practical  method  for  selecting  target  elevation 
estimates.  It  also  studies  all  the  candidate  signal  elevation 
estimates  generated  by  the  SVD  method.  The  results  show  that  the 
optical  track  need  not  be  known  in  the  selection  of  the  target 
elevation  estimates,  if  the  target  range  was  less  than  6.5  km.  They 
also  reveal  the  presence  of  many  false  tracks  branching  out  from 
the  target  and  image  tracks.  A  possible  reason,  supported  by 
computer  simulations,  for  these  tracks  is  given. 

The  presentation  is  organized  as  follows.  Initially,  the 
multipath  radar  data  are  briefly  described  in  Section  II.  Next,  the 
target  elevation  estimates  calculated  with  the  SVD  method  are 
discussed  in  Section  III.  The  presence  of  false  tracks  is  noted, 
and  a  technique  to  select  target  elevation  estimates  derived.  Then, 
computer  simulations  are  presented  in  Section  IV  to  demonstrate 
that  false  tracks  branching  out  from  the  target  and  image  tracks 
can  be  generated  by  adding  array  calibration  phase  errors  to 
uncorrupted  simulated  data.  Finally,  some  concluding  remarks  are 
given  in  Section  V. 
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2.0  MULTIPATH  RADAR  DATA 


The  multipath  radar  data  were  measured  by  STL  in  1981  at 
the  Fraser  Gunnery  Range  near  Portsmouth,  United  Kingdom.  Separate 
vertically  polarized  transmit  and  receive  antennas  were  used.  The 
transmitter  operated  at  a  frequency  of  9.6  GHz,  corresponding  to  a 
wavelength  of  3.125  cm.  The  receive  antenna  was  a  vertically 
oriented,  uniform,  linear  array  comprising  eight  sectorial  H-plane 
horns  stacked  on  top  of  each  other.  The  height  of  each  horn  was 
12.50  cm,  thus  giving  an  element  spacing  of  4.0  wavelengths,  an 
antenna  aperture  of  100  cm,  and  a  beamwidth  of  1.79  degrees  in  the 
boresight  direction.  An  optical  tracker  was  mounted  near  the 
receive  antenna  to  provide  an  independent  measure  of  the  target 
elevation. 

We  have  access  to  five  data  files.  The  first  two  were 
measured  on  the  morning  of  18th  June,  1981,  with  a  Hunter  aircraft 
flying  towards  the  receive  antenna  at  a  nominal  altitude  of  46  m. 
The  height  of  the  antenna  centre  was  6.7  m  above  water.  The  other 
three  files  were  measured  on  the  afternoon  of  25th  June  with  a 
Canberra  aircraft  at  an  altitude  of  61  m.  The  corresponding  height 
of  the  antenna  centre  was  7.9  m.  A  calibration  table,  measured  on 
July  2  with  a  near field  source,  was  provided. 

It  was  noted  in  [1]  that  there  were  date-dependent  offsets 
in  the  optical  tracker  output  and  antenna  boresight  directions,  and 
that  these  offsets  would  result  in  radar  elevation  errors  equal  to 
0.141  and  0.169  beamwidths  for  the  data  collected  on  June  18  and 
25,  respectively.  There  were  also  uncertainties  in  the  horizontal 
reference  direction  for  the  optical  tracker,  and  this  tracker  had 
been  readjusted  before  the  data  measurement  on  June  25. 

Further  details  on  the  experiments  appear  in  [l]  and  [4]. 
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3.0  ELEVATION  ESTIMATES  CALCULATED  WITH  THE  STL  DATA 


We  describe  the  calculation  of  target  and  image  elevation 
estimates,  and  then  present  the  results  and  observations. 

3.1  Calculation  of  signal  Elevation  Estimates 

The  calibration  table  is  applied  to  the  array  snapshots  in 
the  multipath  radar  data.  Each  snapshot  is  then  used  in  a  two-step 
procedure  to  calculate  two  signal  elevation  estimates.  Firstly,  the 
root  variant  of  the  SVD  method  is  used  to  calculate  a  set  of  six 
candidate  signal  elevation  estimates  from  the  snapshot,  using  K=6 
for  the  filter  order  and  M=2  for  the  number  of  signals  present. 
Secondly,  the  power  associated  with  each  candidate  signal  elevation 
estimate  is  calculated  and  the  two  strongest  estimates  are 
identified  as  signal  elevation  estimates.  The  details  of  these  two 
steps  are  given  in  the  appendix. 

The  filter  order  is  calculated  with  the  relation  K=3N/4, 
where  N  is  the  number  of  array  elements.  This  relation  is 
recommended  by  Tufts  and  Kumaresan  in  [2].  The  choice  of  M=2  for 
the  number  of  signals  present  follows  [4],  where  the  best  results 
were  obtained  with  this  number. 

3.2  Results  and  Observations 

The  two  signal  elevation  estimates  (because  M*2)  from  each 
array  snapshot  in  the  five  STL  files  are  plotted  as  individual  dots 
in  Figs,  l  to  5.  No  attempts  are  made  to  link  them,  or  to  associate 
a  dot  with  the  aircraft,  the  aircraft's  image,  or  a  false  estimate. 
The  solid  lines  in  the  figures  are  derived  from  the  optical  tracker 
outputs . 


There  are  several  observations.  Firstly,  the  dots  show  a 
distinct  tendency  to  form  patterns  strongly  suggestive  of  an  output 
"target”  track  close  to  the  optical  track  and  an  output  "image" 
track  below  the  horizontal  direction.  Secondly,  many  false  tracks 
branch  out  from  the  target  and  image  tracks.  They  sometimes  link 
the  target  and  image  tracks  together.  Examples  of  linkages  are 
clearly  visible  at  ranges  3.4,  4.2  and  6.5  km  in  Fig.  5.  Thirdly, 
the  output  target  tracks  are  biased.  The  biases,  or  the  distances 
above  the  optical  tracks,  are  about  0.15  and  0.10  beamwidths  for 
the  data  collected  on  June  18  and  June  25,  respectively. 

Fig.  6,  derived  from  the  first  STL  file,  is  typical  of  the 
new  results  obtained  when  only  the  signal  elevation  estimate  with 
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the  higher  elevation  is  plotted.  Below  6.5  km,  the  estimates  form 
an  output  track  close  to  the  optical  track.  Based  on  this 
observation,  a  target  elevation  estimate  can  be  obtained  from  an 
array  snapshot  in  three  steps:  calculate  a  set  of  six  candidate 
elevation  estimates  with  the  snapshot;  find  the  two  with  the 
highest  power;  and  pick  the  one  with  the  higher  elevation.  This 
procedure  is  more  practical  than  the  one  in  [4],  because  an  optical 
track  is  usually  unavailable. 

The  bias  in  output  target  elevation  estimates  confirms  the 
presence  of  offsets  in  the  optical  tracker  output  and  antenna 
boresight  directions.  It  is  date-dependent,  because  the  optical 
tracker  had  been  readjusted  before  STL  files  3  to  5  were  measured 
on  June  25. 


4 


ELEVATION  (BN) 


RANGE  (KM) 


Fig.  2  Signal  elevation  estimates  calculated  with  the  singular 
value  decomposition  method  and  the  second  STL  data  file. 
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Fig.  3  Signal  elevation  estimates  calculated  with  the  singular 
value  decomposition  method  and  the  third  STL  data  file. 
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Fig.  4  Signal  elevation  estimates  derived  from  the  fourth  STL 
data  file. 


8 


ELEVATION  (BW) 


RANGE  (KM) 


Fig.  5  Signal  elevation  estimates  calculated  with  the  singular 
value  decomposition  method  and  the  fifth  STL  data  file. 
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Fig.  6  The  larger  of  the  two  signal  elevation  estimates  derived 
from  the  first  STL  data  file. 
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4.0  EFFECT  OF  CALIBRATION  PHASE  ERRORS  ON  ELEVATION  ESTIMATES 


We  have  studied  the  effect  of  calibration  phase  errors  on 
elevation  estimates.  The  following  results  are  rather  general: 

(a)  The  direction  estimates  arc  usually  biased.  This  bias  is 
dependent  on  the  calibration  phase  errors,  and 

(b)  The  output  target  and  image  tracks  usually  have  kinks  or 
branches.  The  details  are  also  dependent  on  the  calibration 
phase  errors. 

Some  examples  are  given  in  the  remainder  of  this  section  to 
illustrate  the  above  observations  and  to  show  that  calibration 
phase  errors  can  explain  some  of  the  features  in  Figs,  l  to  5.  The 
antenna  and  target  geometry  in  the  computer  simulations  are  the 
same  as  those  of  the  STL  experiment  on  June  18.  The  reflection 
coefficient  of  the  water  surface  is  equated  to  -0.794, 
corresponding  to  2.0  dB  in  power  attenuation  and  180  degrees  in 
phase  change.  Target  signal-to-noise  ratio  (SNR)  at  the  array 
elements  is  modelled  to  be  inversely  proportional  to  the  fourth 
power  of  range,  so  that  a  50  percent  reduction  in  range  leads  to  a 
12  dB  increase  in  SNR.  A  set  of  321  array  snapshots  are  generated 
as  the  target  range  decreases  from  10  to  2.0  km.  Calibration  phase 
errors,  if  needed,  are  those  given  in  the  following  table. 
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ELEMENT 

PHASE  ERROR  (DEG) 

1 

4,769692 

2 

-2.650831 

3 

1  -9.087350 

4 

-6.525987 

5 

-5.240850 

6 

10.306355 

7 

8.923951 

8 

-0.494981 

Table  of  calibration  phase  errors  applied  to 
simulated  data. 


They  are  defined  as 

phase  with  error  =  phase  without  error 

+  calibration  phase  error  . 

The  average  error  magnitude  in  the  table  is  6.0  degrees,  and  the 
largest  is  about  10.3  degrees.  In  the  convention  used,  the  first 
array  element  is  the  bottom  element  in  the  array. 

Fig.  7  serves  as  a  reference  for  comparison.  It  is  obtained 
by  plotting  the  six  candidate  signal  elevation  estimates  (because 
the  filter  order  is  given  by  K=6)  calculated  with  each  array 
snapshot  and  the  SVD  method.  Calibration  phase  errors  are  not 
applied  to  the  data,  and  the  SNR  is  26  dB  when  the  target  range  is 
10  km.  The  estimates  produce  six  tracks:  a  target  track  close  to 
and  above  the  zero  beamwidth  elevation,  an  image  track  close  to  and 
below  the  zero  beamwidth  elevation,  and  four  false  tracks  with 
approximately  constant  elevations  at  -3.0,  -1.8,  1.8,  and  3.0 
beamwidths . 

The  false  tracks  in  the  above  figure  can  be  removed  by 
plotting  only  the  two  strongest  estimates  obtained  from  each  array 
snapshot.  Not  visually  observable  in  the  figure  is  the  bias  in 
direction  estimates.  This  bias  makes  the  target  and  image  elevation 
estimates  larger  than  the  correct  values  by  an  average  of  about 
0.038  degree,  or  0.02  beamwidth. 

Fig.  8  is  obtained  by  applying  the  calibration  phase  errors 
in  Table  1  to  the  data.  As  the  target  range  decreases,  the  lowest 
false  track  moves  up  to  the  position  of  the  next  lowest  false 
track,  which  in  turn  moves  up  to  the  position  of  the  image  track, 
and  so  on.  This  upward  movement  eventually  stops,  and  kinks  appear 
in  the  output  target  and  image  tracks.  The  amplitudes  of  the  kinks 
decrease  with  the  target  range. 

Fig.  9  is  obtained  by  reversing  the  signs  of  the 
calibration  phase  errors.  The  tracks  move  downwards  as  target  range 
decreases . 

Fig.  10  is  obtained  by  reducing  the  SNR  to  6.0  dB  at  10  km 
target  range  and  plotting  only  the  two  strongest  elevation 
estimates  obtained  with  each  array  snapshot.  This  figure  has  many 
of  the  kinks  and  track-branching  features  found  in  Figs.  1  to  5. 
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Fig.  7  Candidate  signal  elevation  estimates  derived  from 

simulated  data.  SNR  for  10  km  target  range  is  20  dB. 
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Fig.  8  Candidate  signal  elevation  estimates  derived  from  the 
simulated  data  used  in  Fig.  7.  Table  1  for  calibration  phase  errors 
applied. 
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Fig.  9  Candidate  signal  elevation  estimates  derived  from  the 
simulated  data  used  in  Fig.  7.  The  calibration  phase  errors  are  the 
negative  of  those  in  Table  1. 
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Fig.  10  Signal  elevation  estimates  derived  from  simulated  data. 
SNR  for  10  km  target  range  is  6  dB.  Table  1  has  been  applied  to  the 
data . 
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5 . 0  CONCLUSIONS 


We  have  carried  out  a  new  analysis  of  the  multipath  radar 
data  measured  by  STL  and  observed  two  new  results.  Firstly,  we  can 
get  a  target  elevation  estimate  from  an  array  snapshot  without 
using  the  optical  tracker  output,  if  the  target  range  was  less  than 
6.5  km.  The  procedure  is  to  calculate  a  set  of  six  candidate  signal 
elevation  estimates  with  the  SVD  method,  find  the  two  with  the 
largest  power  estimate,  and  then  pick  the  one  with  a  higher 
elevation.  Secondly,  many  false  tracks  branch  out  from  the  target 
and  image  tracks.  Tracks  with  similar  properties  can  be  generated 
by  adding  antenna  calibration  phase  errors  to  uncorrupted  simulated 
data. 


The  presence  of  false  tracks  branching  out  from  the  target 
and  image  tracks  indicates  that  accurate  measurement  of  the 
calibration  phases  of  a  low-angle  tracking  array  antenna  is  crucial 
when  the  SVD  method  is  used  and  the  array  antenna  has  a  small 
number  of  elements.  Otherwise,  long-range  data  may  produce 
misleading  results  on  target  manoeuvres.  For  example,  a 
horizontally  flying  target  may  be  wrongly  interpreted  as  a  target 
flying  upwards. 
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APPENDIX 


CALCULATION  OF  ELEVATION  ESTIMATES 


Presented  are  the  calculations  of  candidate  signal 
elevation  estimates  and  signal  elevation  estimates.  There  are  three 
sections.  The  first  contains  the  notations  and  signal  model  used; 
the  second,  the  method  for  calculating  a  set  of  K  candidate 
elevation  estimates  from  an  array  snapshot;  and  the  third,  the 
culling  of  M  signal  elevation  estimates  from  the  set  of  K  candidate 
signal  elevation  estimates. 

A.l  Notations  and  Signal  Model 


The  array  elements  are  labelled  n=l  to  N,  with  the  first 
element  at  the  bottom.  An  array  snapshot  is  denoted  by  vector  y, 
where  y=(y^,y2»  •  •  • aud  y^  is  the  output  of  the  n-th  array 
element.  Superscript  T  denotes  the  transpose.  An  array  steering 
vector  for  elevation  <f>  is  denoted  by 


v(«)  =  z 


(1/2, Z^ 


.  ,Z 


M'1 


(A.l) 


Where 


z  =  expt  j2rrd  sin  <#>]  ,  (A. 2) 

and  d  is  the  array  element  spacing  in  wavelengths. 

The  signal  model  is 

y  =  V  a  +  ,  (A. 3) 


where 


(V(0,)  ,v(02)  , . . .  ,v(0^)  (A. 4) 

is  an  NxK  matrix  of  array  steering  vectors,  K  is  the  number  of 
signals,  {0, ,02' •  •  are  the  signal  directions, 

a  =  (a^ , a2, . .  •  / a|^)  (A.  5) 

is  a  vector  for  the  complex  signal  amplitudes  at  the  array  centre, 
and  If  is  the  noise  amplitude. 

Estimates  of  the  complex  signal  amplitudes  are  calculated 
with  (A. 3),  which  gives 
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where  superscripts  H  and  #  denote  the  Hermitian-transpose  and  the 
Moore-Penrose  pseudo inverse,  respectively. 

A. 2  Computation  of  candidate  Elevation  Estimates 

There  are  three  steps  in  the  calculation  of  candidate 
elevation  estimates. 

Step  1 , 

Calculate  a  set  of  K  filter  coefficients  with  the  SVD  method. 
Denote  them  by  • 

Step  2. 

Construct  a  polynomial  as 

F(2)  =  1  +  g^2-’  +  +  ...  +  .  (A. 7) 

Step  3. 

Find  the  K  complex  roots  of  F(z)  ,  identify  them  as 
{z. ,  Zj, . . . ,  Zg} ,  and  calculate  the  candidate  elevation  estimates 
with  the  relation 

arg(z,^)  =  2wd  sin  0^  ,  k=l,2,...,K  ,  (A.  8) 

where  z^,  is  the  k-th  root  of  F(z)  and  arg(2)  is  the  phase  of  z^. 
in  radians. 

A. 3  Culling  of  Signal  Elevation  Estimates 

The  M  signal  elevation  estimates  are  defined  as  the  M 
candidate  signal  elevation  estimates  with  the  largest  power 
estimates.  The  power  estimates  are  calculated  with  (A. 6)  in  three 
steps . 

Step  1. 

Find  the  eigenvalues  and  eigenvector  of  the  K  by  K  matrix 
product  v'V,  Arrange  them  in  the  order  of  decreasing 
eigenvalues.  Denote  them  by  {^-,^2, . . .  and  {u^u,, .  . .  ,u„} , 

respectively. 

Step  2 . 

Calculate  the  pseudoinverse  of  v”v  as 


(A. 9) 


(y«  Y)#  ^  ^-1 

lt*1 

The  prime  in  the  summation  indicates  that  only  the  terms  with 
eigenvalues  larger  than  are  included. 

Step  3 

(a)  Calculate  vector  ft  with  (A. 6). 

(b)  For  k=l  to  K,  calculate  the  power  associate  with  candidate 

elevation  estimate  <p^  from  the  components  of  ft  as  |  |  ^ . 

(c)  Identify  the  M  canaidate  elevation  estimates  with  the 
largest  power  as  the  desired  M  signal  elevation  estimates. 

The  choice  of  the  threshold  value,  i.e.  in  Step  2  is 

based  on  the  results  of  a  separate  study  which  indicate  that  the 
best  target  and  image  tracks  are  obtained  if  this  value  is  used  to 
process  the  STL 
data. 
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